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Abstract. The parabolic Af-membranes problem for the p-Laplacian and the 
complete order constraint on the components of the solution is studied in what 
concerns the approximation, the regularity and the stability of the variational 
solutions. We extend to the evolutionary case the characterization of the 
Lagrange multipliers associated with the ordering constraint in terms of the 
characteristic functions of the coincidence sets. We give continuous dependence 
results, and study the asymptotic behavior as t — > oo of the solution and the 
coincidence sets, showing that they converge to their stationary counterparts. 



Dedicated to V.A. Solonnikov, on the occasion of his 75th birthday, with admi- 
ration and friendship. 



1. Introduction 

The aim of this work is to analyze the quasihnear parabolic A^-system associated 
with the scalar operator involving the p-Laplacian in the elliptic part 

(1.1) Pu, = dtu, - V ■ {\Vu,\P-^\/u,) , i = l,...,N, 

with 1 < p < oo, dt — d/dt and V — {d/dxi, . . . ,d/dxd), in a space-time cylin- 
der fix = X (0,T), n C K'^, in the case in which the solution u = u{x,t) = 
(ui, . . . , un) has all its components completely ordered 

(1.2) ui > U2 > ■ ■ ■ > UN , a.e. {x,t) G fir, 

and subjected to a given nonhomogeneous term / = f{x,t) — (/i, . . . , f^) and 
given boundary conditions. For simplicity, we assume 

(1.3) M = on Et 9fi X (0, T) and u^h on = O x {0}, 

for given Cauchy data h. 

The time independent case corresponds to the classical A'^-membranes problem 
which can be formulated as an elliptic variational inequality. It has been studied 
for different types of operators (see [20l [211 HI [2 [3] ) associated with a convex subset 
of a Sobolev space determined by the constraint p.2p . In the recent papers [H [3] it 
has been shown, in particular, that the A^-membranes problem can be interpreted 
as a reaction-diffusion system with additional discontinuous nonlinearities. In the 
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evolutionary case (II. ip . it will be shown in this work that the solution u solves a 
parabolic system of the form 

(1.4) Pu ^ f + R{x,t,u) in f^T, 

where Pu = {Pui, . . . , Pun) and each of the components of the nonlinear reaction 
term R depends on (x, t) through linear combinations of the /j, 1 < i < iV, and on u 
through the characteristic functions ^j,k = ^j^k{x, t) of the N{N — l)/2 coincidence 
sets 

(1.5) Ij,k^{{x,t)enT : Uj{x,t) = ...^Ukix,t)} , l<j<k<N, 

i.e., ^j,k{x,t) = 1 if (.T,<) G Ij^k and Xjj^[x,t) = otherwise. 

We can illustrate the general form of the system (|1.4p for iV = 3 (see [2]) 
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= h + 




-/l)^l,2 
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-h 


-/l)^l,3 
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-,/2 
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which contains the simpler case N = 2, that corresponds to the first two equations 
with ^2,3 = and ^1^3 = 0, in which case the third equation is independent of the 
first two. Noting that, in general, ^j_k — ^j.j+i^j+i,j+2 ■ ■ ■ ^k-i,k, for 1 < J < fc < 
N , in (II. 6p the last terms containing ^1^3 — ^1^2^2,3, are in fact doubly nonlinear 
in u. This introduces additional difficulties in analyzing the stability of the system 
with respect to the perturbation of the data. In fact, in section 3, we show that 
the sufficient conditions on the averages of the components of /, obtained in [3] for 
the stability of the coincidence sets Ij^k in the stationary problem, extend to the 
parabolic case as well. In particular, for N = 'i, they take the form 

/i ^ /2, /2 ^ h, fi + /3 ^ a.e. in Qt- 

We notice that the stability result on the X^.j. is not a direct consequence of the 
stability of the solution u with respect to the data / and h, which can, however, 
be obtained by direct variational methods, as we also show in subsection 12.41 

Classical monotonicity methods (see |T3] , for example) or the theory of accretive 
operators and evolution inclusions in Banach spaces (see [T2], [TH], [H] or [T] and 
their bibliography) are directly applicable and yield general results on the exis- 
tence of solutions to our problem, when formulated as a variational inequality in 
the convex set associated with the constraints (|1.2p . In section 2, we introduce 
an approximation of the variational inequality formulation and we obtain directly 
useful a priori estimates for the existence of solutions. We remark that we as- 
sume the p'-integrability of / and rely on the p-integrability of a compatible h and 
its derivatives, but we do not require the boundedness of h nor of the variational 
solution globally in fix- 
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Considering the relation of the upper and lower membranes (in particular, the 
two-membrane problem) with the obstacle problem and of the inner membranes of 
the A^-problem, with > 3, with the two-obstacles problem, we apply the dual 
estimates for unilateral parabolic problems (see [6], [12] or [11]) to obtain Lewy- 
-Stampacchia type inequalities 

i N 

(1.7) /\ /, < Pu, < V fj a.e. in r!T, i = l,...,N, 

for the parabolic operator (jl.ip . Here we use the notation 

fc fc 

=CiV...Vefe =sup{6,...,6} and /\ = 6 A . . . A^- = mfUi, • • • , 6} 

i=l 1=1 

and we also denote ^+ = ^ V and = — (^ A 0). 

We also show how the estimates on Pui imply that the variational solution to 
the A^-membranes problem solves a.e. a system of the type (|1.4I) . for an explicit 
R with the same p'-integrability as /, extending the analogous result obtained in 
[3] for the stationary problem. This implies, in particular when / is bounded, the 
Holder continuity of the solution and of its gradient. In fact, this is an immediate 
consequence of known estimates for the parabolic operator (jl.ip . even without 
knowing the explicit form of R, as we observe in section 2. Even for the linear case 
p — 2, for which we can apply Solonnikov's estimates in W^'^(S7t), the regularity 
obtained here for the solution of the evolutionary A'^-membranes problem is new. 

In section 3, we study the asymptotic convergence, when t — > oo, of the solution 
u{t) to the corresponding solution of the stationary problem of [3 in i^(ilT) (here 
we denote Z,^(riT) = [-^^(^^t)] ), in the case p > 2. We show how a modest 
convergence of the solution, obtained as in [19] . also implies the asymptotic sta- 
bilization of the evolution coincidence sets towards the stationary ones, under a 
natural nondegeneracy assumption identified in [3]. 

Finally, we observe that most results still hold, with suitable adaptations, for 
more general quasilinear parabolic scalar operators 

Pu ~ dfU — V • {a{x, t, Vu)) , 

in particular, for strongly monotone vector fields a(-,^), with p-structure as in [3], 
as well as more general data / in L'^{Q,T; L'^(ri)) (see [4]). 

For simplicity of presentation, we limit ourselves here to the case of the p- 
-Laplacian with homogeneous Dirichlet data, i.e., we consider only variational so- 

1 r 1 1^ 

lutions in the usual Sobolev space Wg^^iQ) = IWo'^in)] , for 1 < p < oo. The 
case of a time-dependent Dirichlet boundary condition is more delicate and will be 
considered in [ITl. 
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2. Approximation and regularity of variational solutions 

Let i7 C M'' be a bounded domain with a Lipschitz boundary, let T > 0, and 
define the space-time domain Jly := fl x (0,T), with parabohc boundary dpflx ■— 
St U Hq. We use A^- vectorial notation for vector fields 

w := (wi, . . .,wn) e 

and function spaces F := [F]^ . For 1 < p < oo, define the differential operator 

VpW — {VpWi, . . . , VpWN) , VpWi :— |Vwi|^ ^ Vwi and ApW = V • VpW. 

We assume the data satisfy 

(2.1) feL^'iVLr) and /ieKni^(f7), 

where p' ^ p/[p — 1) and IK is the closed convex subset of W\'^\Vl) defined by 

(2.2) IK = e Wl'^{n) : vi>...>VN, a.e. in f^j . 

2.1. Variational formulation of the problem. The evolutive iV-membranes 
problem for the p-Laplace operator consists in finding a vector field u = u{x,t) 
such that 

(2.3) w e LP (O, T; W^^m) n C ([0, T]-L\n)) , 

(2.4) dtueLP'{Q,T;W-^^P\n)), 

(2.5) u{t) e K, a.e. t e (0, T), m(0) ^he L\VL), 
and, for a.e. t G (0, T) and all t> G IK, 

(2.6) (dtu{t),v - u{t)) + [ Vpuit) : V (-y - M(i)) > / /(<) ■ (v - uit)) . 

Here, (•, •) denotes the sum of the N duality pairings in W'^''^ (Vl) x Wl'^{n) of 
the components of the vector fields, and A : B denotes the scalar product of the 
matrices A and B. 

We observe that, by a simple comparison argument, there exists at most one 
solution of (|2.3p - (|2.6p . the variational inequality formulation of the evolutionary 
A^-membranes problem. 

2.2. The approximating problem. We approximate the variational inequality 
(|2.6p using a bounded penalization. For that purpose, for each e > 0, let 6'e be the 

e<-e 

-e <0 <0 

e>Q. 



real function defined in [— oo, +oo] by 

e,{e) = <^ e/e if 

if 
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The approximating penalized problem is the system of boundary value problems 
defined as follows: 



(2.7) 



uf = on YiT and ttf = hi on VLq 



with i = 1,...,7V, and the convention uq = +oo and wat+i = — oo, where for 
i^l,...,N, 

(2.8) eo-niax{^i±— ±^ :* = 1,...,7V}, ^ ^ - (/i + ■ • • + /.), 
(see [3]). Notice that, for i = 1, . . . , A^, we have > and e (il). 

Lemma 2.1. Using the convention vq = +oo and vn+i = —oo, the operator 
N r 

{Bv,w) iiiOeiVi - Vi+i) - ^i-i0e{Vi-i - Vi))Wi, V,weWl'^{n), 

is T-monotone, i.e., 

{Bv-Bw,{v-w)+)>0, yv,weWl'''{n). 
Proof. Since we can rewrite 

N-l „ 

{Bv,w) = V / CiOe - Vi+i) {w., - Wi+i), 

it is enough to observe that 

{Bv — Bw, {v — w)^) = 

N-l „ 

/ ^i{Oeivi - Vi+i) - 6eiw^ - Wi+i)j ({v, - Wi)+ - {v,+i - Wi+i) + ). 
1=1 -^^^ 

As > 0, for i — 1,...,N and Og is monotone nondecreasing, the conclusion 
follows. □ 

Proposition 2.2. Under assumption (|2.ip . the approximating problem p.7p has a 
unique solution (uf,. . .,u%)e 2.^(0, T; T^o'^(f])) C {[0,T]; L'^{n)) such that 

(2.9) i^2,...,N. 

Proof. The existence and uniqueness follow, respectively, from standard results 
concerning monotone operators and comparison (see [13] or [H]), for instance, 
using the Facdo-Galcrkin approximation. We notice that, since f € (^t) C 
LP'(0,T; VF"i'P'(r2)), we obtain, in particular, that dtW £ LP {0,T;W~'^-p' (n)). 
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To prove inequality p.9|) . multiply both the i-th and the (i — l)-th equations by 
{uf — uf_^ — subtract and integrate over Q, obtaining 

I l«- / (Vp<-Vp<_i) •V«-uti-e)+ 

= / [{h~iM<-uU^)+£.^^le,{uU-ul-e)){ul-uU)^ 
Jn 

< / ((/.-/.-i) + (6-e«-i)-(6-i-e.-2))«"<-i-^)+<o. 

Integrating between and t, using the fact that hi > ■ ■ ■ > hpj and the inequality 
^ (Vp< - Vp<_i) • V« - uti - e)+ > 0, 



we get 
(2.10) 



{um-<-i(t)-ef 



<0, 



and so uf < uf_i + e a.e. in VLt- 



□ 



2.3. Existence of variational solutions. The proof of the existence of solution 
for the variational inequality (|2.6p will be done passing to the limit in e ^ 
on the sequence of approximating solutions it^, by using the following a priori 
estimates that can be rigorously obtained through the respective Faedo-Galerkin 
approximations. 

Proposition 2.3. Under assumption (j2.ip . the solution of the approximating prob- 
lem (|2.7p satisfies the following estimates, for a nonnegative constant C , indepen- 
dent of e: 

(2.11) II<IIl-(o,T;L=(o)) + I|V<|Up(o,) < c, 

(2-12) ll^t"f IIlp'(o,T;1v-1'J''(o)) - 

/i < Pu\ < /i+Ci 

(2.13) I f^-i.-l < Puf < /, + ^, (i = 2,...,7V-l) 

, fN-^N-i < Pu% < fN a.e. infix. 
Proof. For each i = 1, . . . , iV, we easily conclude (j2.13p from (|2.7p in the form 

P<=f^+gl m^T, 



where 
(2.14) 



is uniformly bounded in (^t)- 
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Then, multiplying each equation in (|2.7|) by uf and integrating on fit — fix (0, i), 
we get 

Using Poincare inequality, we find 



(2.15) / \um\'+ / |V<r<Co, 

Jn Jut 

where the constant Cq only depends on 11/1111,2(0) and ||/||£,p'(Ot)' H^^^^^i from 
(|2.15p . we immediately obtain (|2.1ip . So 

(2.16) Ap< is bounded in L^' {0,T]W-^'P\n)) independently of e, 

and we conclude (|2.12p by recalling (|2.13p . □ 

Theorem 2.4. Under assumption (j2.ip . the problem (|2.5p - (|2.6p has a unique vari- 
ational solution u in the class (|2.3p - (|2.4p . 

In addition, u'^ ^ u strongly in LP(0, T; VFg'''(ri)) and 

(2.17) Pu" Pu in V^'i^r) - weak. 

Proof. If {u^}e is a sequence of solutions of the approximating problems (|2.7p . by 
the a priori estimates (|2.1ip and (|2.12p . we can extract a subsequence such that, 
as e — > 0, 

in LP(0,T; Wl'^{Vl)) - weak, 
dtu'^dtu in LP'(0,r;Ty-i'?''(f7)) -weak, 

and, by compactness, also u'^ u strongly in U'{U,t). 

Let V e LP{0,T;Wl-P{n)) be such that dtv e Lp' {Q,T;W-^'P' {n)), v{t) e K, 
for a.e. t £ (0,T), and v{0) = h. As {Bv{t),v{t) - u{t)) = 0, we have 

(9t^t^ i; - M^) + /" VpW^ ■.V{v^u') > [ f ■ {v - u") 

It follows from integration by parts that 

{dtu',v~u')^{dtv,v~u')~\ [ \u'{T)-v{T)\^ 

and, using the monotonicity, we get 

{dtv,v-u')+ f \7pV:V{v-u')> f f ■ [v - u') + \ ( \u' {T) ^ v{T)\'' 
Jnr Jnr ^ Jn 



Jnr 

Letting e — > 0, we obtain 

(2.18) {dtv,v-u)+ VpV:V{v-u)> f-{v-u). 

Jnr Jnr 

Now, let w = u + 9{v — u), 6 E (0, 1]. The verification that w can be chosen as 

test function in (|2.18p is immediate. So, 



-/ Vp{u + 0{v~u)) ■.0V{v-u)> ef-{v-u) 
Jqt Jnr 
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Dividing both members by 9 and letting 6* — + 0, we see that u solves the problem 
(9tM, V -u) + / VpM : V(v ~u) > f ■ {v -u), 

for all V such that v e LP{0, T; VTq ^(17)), v{t) e K for a.e. t £ (0, T) and v(0) = /i. 
Using standard arguments (see [E]), also 

(dMt), v{t) - u{t)) + f Vpu{t) : V{v{t) - uit)) > f fit) ■ iv{t) - u{t)), 
Jn Jn 

for a.e. t e (0, T), for all v such that v e LP{0, T; VFq ^(^2)) and e K. 

In order to conclude ()2.17p it is sufficient to recall the estimates ()2.13p for Puf 
and that Vpwf —>■ S7pUi in an appropriate sense. In fact, recalling (j2.14p and using 
equation (|2.7p . we conclude that 



< limsup 



limsup / VpU^ ■ V(tt'^ — u) 

I {f + g')-{u'^u)-{dtu',u^u') 
Jn-r 



0. 



/SIt 

By well-known results (see, for instance, [5 ) this is sufficient to show that — > it 
strongly in LP{0,T; WQ'^{fl)) (notice that ^ g weakly in [VLt), for some g). 

□ 

Remark 2.5. If we assume also that f G i^(f7y), which is a consequence of (|2.ip 
if I < p < 2, the Faedo-Galerkin approach yields directly the regularity 

(2.19) u £ H\0, T; L^inr)) n L°°{0, T; Wl'P{n)) 

through multiplication of (|2.7p by dtuf. 



2.4. Strong continuous dependence. 

Theorem 2.6. Let u* be the variational solution to (|2.5p - (|2.6p corresponding to 
data f* and h* satisfying also (|2.ip and denote 

with q — p'Kl (i.e. q — p' if p > 2 and q ^ 2 if p < 2). Then there exists a positive 
constant c = c{T,p) such that 

(2.20) sup [ \u*{t)-u{t)f + [ |V(m* -m)|^' < ce* ifp>2, 

(2.21) sup j \u*{t)-u{t)\^ + ( j \V{u* -u)\Py <ce* ifl<p<2. 
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Proof. Let V = u*{t) in ()2.6p with data / and h, and v = u{t) in (|2.6p with data 
/* and h* . In the latter case, we have 



(2.22) {dtu*{t),u{t)-u*{t))+ I VpU*{t) -.Viuit) -u*{t)) 

Jn 

> f r{t)-{u{t)-u*{t)). 

From (|2.6p and (|2.22p . integrating between and t, we obtain 

(2.23) i / \u*{t) -u{t)\^ + f f (VpM* - VpM) : V(m* - m) 



1 



< / (/*-/).(„*-„) + - / \h*-h\^ 



In the case p > 2, since 



f f {VpU* -Vpu) -.Viu* -u)>Cp f f\V{u*~u)\P, 
Jo Jn Jo Jn 

the conclusion follows easily by using Holder and Poincare inequalities. 

In the case 1 < p < 2, from (I2.23P we find 

/ \u*{t)-u{t)\^ <e* + f I 
Jo Jo JO 

which, by Gronwall inequality yields, first 

sup / \u*{t) -M(t)P < 
<t<T Jn 



o< 

and, afterwards 



(2.24) / /" (VpM* - Vpw) : V(m* -m) < i(l + re^)e*. 

Jo Jo 2 

Next we consider the following reverse Holder inequality: given < r < I and 

r' = if F e U{Vl), FG e L^{n) and ( \G{:x)f dx < oo in VLt, one has 
' Jo 

\F{x)\''dx^' < (^j \F{x)G{x)\d:^ (^j \G{x)f dx 

Letting r — ^ and, for i ~ 1, . . . , N, F ^ I'^i'^* ^ ""i)!^ we get 



(VpU* - VpUi) • V(m* - Ui) > 



> I / |V 

'o; 



Aj (|V<| + |Vu.|)2-p 

K-^Or^ (l^ (|v<| + |Vu,|)^ 

where f2j = {{x,t) G l^t : |Vw*| + \ Vu,\ > 0}. Thus, if we denote 

2-p 
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by 12231), we conclude ([2?2T|) from 

□ 

2.5. Holder continuity and further regularity of the solution. The regular- 
ity of the variational solutions of the evolution TV- membranes problem does not, 
in general, yield their boundedness for 1 < p < d; but, by Sobolev imbedding, the 
solutions are bounded for p > d and even Holder continuous in the space variables 
for eacht £ (0,T). 

However, estimates (I2.13P and (|2.17p imply that, in fact, Pu has the same regu- 
larity in Ht as the data /. Then, if / G L°°{nT), local and global Holder estimates 
for the evolution p-Laplace equation may be directly applied to bounded solutions 
of the iV- membranes problem (see [5], [2] or [lUj)- In order to illustrate these re- 
sults, we assume in addition that h G L°°{^), which also implies that u G L°°{^lT), 
and consequently that u and Vm are locally Holder continuous. Referring to [7] 
and [14j for the boundary and initial regularity in the space of Holder continuous 
functions C", < a < 1, with the standard parabolic norms, we may state the 
following result: 

Theorem 2.7. Suppose f G L°°{Q) and the initial data h £ C (Ti) D K, < a < 
1. Then the solution u G C" (P't), < a' < a < 1, and Vm G C(^^(Ot), for some 
< /3 < 1. If, m addition, 30 G C^''' and SIh £ C^^ (p), Q < (3 < I, then also 
Vu G C^'iyiT), for some < ^' < /3 < 1. 

In case of a linear operator (p = 2) we can apply directly Solonnikov's parabolic 
estimates (see [13], Thm. 9.1 of page 341). 

Theorem 2.8. Let p = 2. Then, for any f G L'^{Q,t), q > 2, the solution u to 
(|2.5p - (|2.6p satisfies u G W^q'ioc(^7")' which implies, by Sobolev imbeddings, that u 
and Vm are locally Holder continuous, respectively for q > and q > d + 2. If, 

2 2 

in addition, ft. G K H W ^9'^(f2), those results can be extended up to the boundary 
dQ G and up to t = 0, i.e., u G Wq^^ijlr) and u, S/u are Holder continuous 
on fix. 



3. The A^-system and its stability 

The A''-membranes problem can, a posteriori, be regarded as a lower obstacle 
problem for ui, a double obstacle problem for Uj,j — 2<j<N — l, and an upper 
obstacle problem for un- This fact has interesting consequences and, similarly to 
the theory of the obstacle problem that we recall briefly for completeness, allows 
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US to characterize the TV-membranes problem as a nonUnear paraboUc system with 
known discontinuous nonhnearities on the right hand side as in (jl.4p . 

3.1. Dual estimates for obstacle type problems. We consider the scalar two- 
-obstacles problem for the nonlinear operator P defined in (jl.ip . with compatible 
Cauchy-Dirichlet data on dpVLx- Let 

(3.1) (^GLf'(f]T), yye <'P(r!))nL2(r!), 

(3.2) ^1, ^2 e iP(0,T;W^i'P(r!)), i/'i > V2 in f^T, V^i > > t/^s on Et, 
and, for j — 1, 2, 

(3.3) dt^, e LP'(0,T;M^-i'P'(r!)), Pi,, e L^' {nr), Vi(0) > > ^2(0) on n^. 

Using the Lipschitz continuous function 6^ defined in subsection 12.21 for each 
e > 0, we may easily show that the problem 

(3.4) Pw^ + Q2ee{w' -i^2)-Ci0e{'4'i-w')^ip in J^t, 

(3.5) ~ Q on and ^ rj on flo, 

where Ci = {Pipi — and ^2 ~ (-P'02 ~ has a unique solution G 

LP{0,T;W^'P{n)) n C{[0,T];L^{n)), with Pw^ e LP^Qt), uniformly in e < 1. 
Similarly to Proposition [221 it is easy to show that 

-02 — £ < < tpi + e a.e. in fix, 

and, when e ^ 0, as in Theorem [131 that w" w strongly in ^^(0, T; Wj'P(r2)), 
where w is the unique solution of the double obstacle problem 

(3.6) weK'^l^ive LP{0, T- W^^^m : Vi > > V2 in ^t}, 

(3.7) f iPw-ip)iv-w)>0, V-yeK^^, a.e. ie(0,T), 

such that w{0) = 77 on fi. The solution w satisfies also 

w e LP{0,T;Wa'P{n))nC{[Q,T];L^{n)) and Pw e Lp' (nr) 

and, arguing as in Proposition 4.1 of [16j . we can state the following important 
property. 

Proposition 3.1. The solution w to (|3.6p - p.7p . under assumptions (|3.ip - p.3p . 

satisfies the parabolic nonlinear equation 

(3.8) Pw = if+ {Ptp2 - v)^X{w=i,2} - (^^1 - f rX{w=^i} a.e. in f^r- 
In addition, we have the Lewy-Stampacchia inequalities 

(3.9) (p~{Ptpi-(py = ipAPipi <Pw< ipyPi)2 = ip+{Pijj2-v)'^ a.e. in Qt 
and the a.e. in fix necessary conditions for contact with the obstacles 

(3.10) {w ^ ipi} C {Pi^i < Lp} and = V'2} C {PV2 > 
being the inclusions valid up to subsets of Qt with zero measure. 
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Remark 3.2. We note that for the case of only one-obstacle, we have similar 
properties. In fact, if we formally take tjji = +00, we have a lower obstacle problem 

(3.11) w > ^p2 and if < Pw < If \/ P'ijj2 ci.e. in fix, 
and, with ^2 = —00, an upper obstacle problem 

(3.12) w < tpi and If A Ptpi < Pw < ip a.e. in fix- 

Analogously, the semilinear equation holds in each case with the corresponding 
characteristic function, respectively. 

Remark 3.3. As observed in |16j . we have 

(3.13) Pw ^ ip a.e. in {V'2 < w < ipi} 

and due to the fact that both Pw and Pijjj are integrable, we have 

(3.14) Pw — Pijjj a.e. in {w = ipj} for j — 1,2. 

Using the regularity of Theorem 12.41 we easily see that each component Ui of 
the iV-membranes problem solves an obstacle type problem (I3.6l) - (|3.7p with ip = fi, 
■01 = and "02 — Ui^i (with the conventions uq = +00 and u^+i = —00 

corresponding to the one-obstacle problems). Hence, we have from p. lip . p.9p 
and (|3.12p . respectively, a.e. in 17t, 

/l < PUI < flVPU2 

/, APu,_i < Pu, < /,VPu,+i (i = 2,...,7V-l) 

fN /\ Pujsi-i < Pun < /at a.e. in Ot- 
By simple iteration, we have shown the following Lewy-Stampacchia type in- 
equalities, that extend Theorem 3.5 of [3] to the evolution A'^- membranes problem 

Theorem 3.4. The solution u of (I2.5p - (|2.6p satisfies 

i N 

A fj <Pu^<\/ fj a.e. m^T, i^l,...,N. 

j=l j=i 

3.2. The nonlinear iV-system. As a consequence of the equivalence of the A^- 
-membranes inequality with two one-obstacle problems and N ~ 2 two-obstacles 
problems, we may prove the equivalence of this inequality with a A^-system of 
equations, strongly coupled by the ^'-^""^^ coincidence sets Ij^k defined in (II. Sp . 
Indeed, we can argue as in section 4 of [3], and since we know that Pui G LP (Ot), 
for alH = 1, . . . , A^, we have on each coincidence set 

Puj = • • • = Puk a.e. in Ij^k = = = ""*:} 

and we conclude, for each j < i <k, 

Pui = {f)],k a.e. in Ij^k, 
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where we introduce the averages of / by 
/f\ ■^ ^ fk 



I < i <k < N. 



k-j + l 

On the other hand, in the complementary sets ^It \ Ij,k, for each i > k > j or 
i < j < k, we have 

Pui = f, a.e. in fiy \ Ij,k, 
and we conclude, as in ^3 , the following explicit form for (|1.4p . 

Theorem 3.5. The variational solution of the N -membranes problem (|2.5p - (|2.6p 

satisfies the system {i — 1, . . . , N) 

(3.15) Pu,^f,+ Yl ^Af]^3,k a.e. mUr, 

1 < j < k < N,j < i < k 

where ~^j^k denotes the characteristic function of each Ij_k and 

(f)],k - {f)]M^i if i^j 



(f)],k - {f)j + l,k 



if i — k 



K (k^ml-j+i) i(f)w-i - Uf, + fk)) if ]<i<k. 
For the particular case TV = 3 (and iV = 2), we can easily deduce (|1.6p from 



3.3. Convergence of coincidence sets. From Theorem 12.61 we know that if for 
sequences 

(3.16) r >f mi«(f^T), g = p'A2, 



(3.17) 



h in L'^iVtr), 



then, the corresponding solutions of (|2.5|) - (|2.6p also converge 

(3.18) u'' >u inC\[0,T]-L'{^l))r\LP{0,T-Wl'P{^l)). 

Consequently, we have 

Apu"^ ApM in LP'{0,T; W^'^^p' (Q)) - weak 

and, by Theorem 13.41 also 

Pu" --^ Pu in L'^{VIt) - weak. 

Since the characteristic functions '^j ^ = ^{u''.=---=ui} satisfy < X^j, < 1 a.e. in 
fir, there are k ^ L°°{Q,t) such that 

^),k -17^ ^j,fc in L°"{^t) ~ weak * . 



Passing to the limit in 



E ^^'[r]^lk. 

1 < j < k < N,j < i < k 
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we obtain, for each i = 1, . . . , N, 

1 < j < k < N,j < i < k 

which compared with (j3.15|) yields 

E K''[fW,,k-^l,) = a.e. inf^T. 

1 < j < k < N,j < i < k 

Arguing exactly as in the proof of Theorem 4.6 of [3], we conclude, under the 
same nondegeneracy assumption for the limit data, namely 

(3.19) ^ a.e. inr!T, for aU i, j, fc e {1, . . . , TV}, with i < j < fc, 

that ^j^k = -"^j k prove the following stability property for the respective coin- 
cidence sets Ijj, — {Uj — ■ ■ ■ ^ u^}. 

Theorem 3.6. Under the convergence assumptions (j3.16p and p.l7p . the char- 
acteristic functions associated with the convergent variational solutions p.lSp also 
converge 

for any 1 < s < oo, all I < j < k < N , provided the nondegeneracy condition 
([XTO)) holds. 



3.4. Asymptotic stabilization as i oo. In this section we assume p > 2 and 
we consider the unique solution u°° to the stationary A^-membranes problem for a 
given /~ G LP'iny. 

(3.20) u^eK: [ VpM°° : V(v - M°°) > [ /°°-(v-m°°), ^veK. 
Jn Jn 



Supposing that the problem (|2.5p - (|2.6p is solvable for all T < oo and that 
f{t) — > /°° in (il) as t — » oo in the sense 

(3.21) r' f \f{t)- r\''' ^0 ast^oo, 
Jt Jn 

by the results of [12], the evolutive solution u{t) is such that 

(3.22) u{t) > in L^{n), 

t—>-QO 

(3.23) / [ |VM(t) - Vm°°|'' — >0 ast^oo. 
Jt Jn 

By the results of ^3], the stationary solution also solves the nonlinear A^-system 

(3.24) -Apur = /r+ E K^'ifn^^ a.e.inn, 

^ < j < k < N, j <i<k 

where = -^{ti°°= - =Mjf } denotes the characteristic function of the limit coinci- 
dence set /~ ={x€n:u°°{x) = --- = u^(x)}. 

Denoting by ^j,k{t) the characteristic functions of Ij,k{t) — {uj{t) — ■ ■ ■ = Uk{t)} 
at time i, we have the following asymptotic convergence result as i — > oo. 
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Theorem 3.7. Under assumption (I3.21|) . the variational solution of the evolution 
N -membranes problem converges to the corresponding stationary solution in the 
sense (|3.22[) and (j3.23p . In addition, the characteristic functions satisfy 

(3.25) ^]At)^^T-k ast^oo inL'^iQ), 
for any 1 < s < oo, for all 1 < j < k < N , provided we assume 

(3.26) {f°^)^J ^ {f°^)j+i,k a.e.mn, for all I < i < j < k < N. 
Proof We rewrite ((3?23l) for w{t) = u{t) - u°° as 

( f \\/w{T)\PdT^ f f \\/w{t + cr)\Pda — >0 asi^oo, 
Jt Jn Jo Jn 

and this convergence can be interpreted as 

(3.27) w,(i) — >{] &st^ooinLP{Q,l-Wl'P{^)), 
where we define E L°°(0, oo; LP{0, 1; Wl'^ifl))) as 

Wi{t)ia, •) = w{t + a, •) e Wl^P{n), a S (0, 1). 

Consequently, from ()3.27|) we have 

ApMj(<) ^ ApU^ as i oo in Lp' {0,1: W^^'P' {Q)) - weak 

and, recaUing the estimates of Theorem 13.41 and the assumption (|3.21[) . we may 
conclude 

{dtUf - ApMjj)(i) ~ApU^ as i ^ oo in Lp' {0, l;LP'{n)) - weak. 

Since u^{t) solves p. 151) . a.e. in S7 and for a.e. t > 0, we can pass to the limit, 
as t ^ oo, in LP (0, 1; LP (fi)). As in the proof of Theorem 13.61 (and Theorem 4.6 
of [3]), we conclude that assumption (|3.26p implies the convergence ^j,fe(t) — > 
as i — > oo, first as functions of L°°(0, 1; -L°°(il)) with the weak-* topology and, 
afterwards, also in the sense of (|3.25[) . Indeed, since they are characteristic functions 
and any subsequence of ^j.k{t) has the same limit ^'j°k^ their weak convergence 
implies the strong convergence in L'^^fl) for all 1 < s < oo. □ 
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